Derealization and wave-packet dynamics in one-dimensional diluted Anderson models 
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We study the nature of one-electron eigen-states in a one-dimensional diluted Anderson model 
where every Anderson impurity is diluted by a periodic function f(l) . Using renormalization group 
and transfer matrix techniques, we provide accurate estimates of the extended states which appear 
in this model, whose number depends on the symmetry of the diluting function /(/). The density 
of states (DOS) for this model is also numerically obtained and its main features are related to the 
symmetries of the diluting function /(?). Further, we show that the emergence of extended states 
promotes a sub-diffusive spread of an initially localized wave-packet. 

PACS numbers: 63.50.+X, 63.22.+m, 62.30.+d 



I. INTRODUCTION 

The nature of electronic states in disordered tight- 
binding models with site-diagonal uncorrelated disorder 
was firstly studied by Anderson [l| when the localization 
of quantum states was discussed in connection with the 
transport properties of a random lattice. One of the most 
remarkable effects of disorder, which has been demon- 
strated by several authors, is the exponential localization 
of all one-electron eigen-states in the one-dimensional 
Anderson model, irrespective to the strength of disor- 
der |2j. The electron localization in random lattices is 
a result of destructive quantum interference due to in- 
coherent backscattering. Therefore, the localization phe- 
nomenon occurs in general model systems involving wave 
propagation in random media. 

Recently, the existence of delocalized states in sev- 
eral variants of the low-dimensional Anderson models 
has been reported EM SBIllSISIISIilllllll 

[H E El E E3 In these works, the presence 
of short or long-range correlations appears as the fun- 
damental mechanism responsible for the emergence of 
extended states. This theoretical prediction about sup- 
pression of localization was recently confirmed experi- 
mentally in doped polyaniline |2flj , semiconductor super- 
lattices with intentional correlated disorder [2l| as well 
as microwave transmission spectra of single-mode waveg- 
uides with inserted correlated scatterers |23- Among 
these models, the diluted Anderson chain has attracted 
a renewed interest 

!2llHlil EES- Hilke [H intro- 
duced an Anderson model with diagonal disorder diluted 
by an underlying periodicity. The model consisted of two 
interpenetrating sub-lattices, one composed of random 
potentials (Anderson lattice) and the other composed of 
non-random segments of constant potentials. Due to the 
periodicity, special resonance energies appear which are 
related to the lattice constant of the non-random lat- 
tice. The number of resonance energies is independent 
of the system size and, therefore, it was conjectured that 
these states shouldn't have any influence on the trans- 
port properties in the infinite size limit. In ref. the 



authors presented a simple model for alloys of compound 
semiconductors by introducing a one-dimensional binary 
random system where impurities are placed in one sub- 
lattice while host atoms lie on the other sublattice. The 
existence of an extended state at the band center was 
demonstrated, both analytic and numerically The di- 
luted Anderson model was recently extended to include 
a general diluting function which defines the on-site en- 
ergies within each non-random segment p^ |. Using a 
block decimation approach, it was demonstrated that this 
model displays a set of extended states, the number of 
which strongly depends on the length of the diluting seg- 
ments and the symmetries of the diluting function. 

In this work we will use the matrix decimation method 
and the transfer matrix technique to provide accurate es- 
timates of the set of extended states in general ID diluted 
Anderson models. The density of states (DOS) will be 
shown to display a set of gaps whose number and width 
depend on the set of impurities used to dilute the disor- 
dered lattice. Further, we will investigate the influence 
of such resonant extended states on the wave-packet dy- 
namics. Our results suggest that, in spite of the num- 
ber of resonance energies being independent of the sys- 
tem size [23| , these extended states induce a sub-diffusive 
spread of an initially localized wave-packet. 



II. MODEL AND FORMALISM 

The standard one-dimensional Anderson model is de- 
scribed by a tridiagonal Hamiltonian 



H 



J2 e M ><j\ +t J2^ xj+M+y ><i-i|] (!) 



where disorder is introduced on the site energies €j which 
are uncorrelated random numbers chosen from a previ- 
ously defined distribution. In our calculations, we will 
use energy units such that the hopping term t = 1 and 
the random site energies will be taken uniformly from the 
interval [—0.5,0.5]. The diluted Anderson model is con- 
structed by introducing a sequence of k new sites between 
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each original pair of neighboring sites. These sequences 
are all identical and the on-site energies within such se- 
quences are given by f(l), I — 1, 2, k. 

To study the properties of the one-electron eigen- 
states, we employ a transfer matrix calculation (TM) in 
order to obtain the Lyapunov exponent defined as the 
inverse of the localization length L c . The Schrocdinger 
equation for the present tight-binding model is: 



-1 + C n +1 = ECn 



(2) 



where <& = X)n c «l n - > ^ s an eigenstate with energy E. 
The above equation can be rewritten as a transfer matrix 
equation 



Cn+l 
C/v 



E-e r 
1 



-1 




, n = 0,l,2,...,iV. 

Based on the asymptotic behavior of the matrix prod- 
uct ri n _i Qni where Q n are 2x2 transfer matrices, the 
Lyapunov exponent 7 can be defined as: 



7 



lim — log ■ 



n:=iQnz(o)\ 

k(o)| ; 



(4) 



where z(0) = is a generic initial condition. In addi- 
tion, the Lyapunov exponent can be obtained using a dec- 
imation renormalization-group (DRG) technique which 
is based on the particular form of the equation of mo- 
tion satisfied by the Green's operator matrix elements 
[G(E)] hJ = {i\l/(E-H)\jM: 



(E 



^)IG(E)} 



n-\-fi,n 







(5) 



-i[G(E)} n 



+fi—l,n 



t °n+ t i.,n+ t i.+ l[G{E)]n+^+l,n 1 



where e° 



t. After crim- 



en anu — o n _ i >n 

inating the matrix elements associated with a given 
site, the remaining set of equations can be expressed 
in the same form as the original one, but with ener- 
gies and hopping amplitudes renormalized. Therefore, 
the decimation-renormalization technique consists in re- 
moving iteratively the sites 1,2, N — lof the system 
obtaining in that way the effective energies of the exter- 
nal sites and the effective hopping between them by using 
the iterative equations 




FIG. 1: (a) Dashed lines represent the Normalized density of 
states (DOS) as a function of energy E obtained using the 
Dean's method for a chain with N = 10 sites of the diluted 
ID Anderson model with dimer segments of identical on-site 
energies to- The DOS displays two pseudo-gaps for eo > 0, 
which are reminiscent of the gaps present in the DOS of the 
corresponding pure models (shown as solid lines), (b) The 
width of the gap (A) for a diluted chain without disorder 
versus eo- The width A displays a linear behavior A oc ej for 
any length of the diluting segment. 



evidence of the localized/delocalized nature of the one- 
electron states. An oscillatory behavior of t^J\ as a 
function of N , signals an extended state. On the con- 
trary, an exponentially localized state results in an expo- 
nential decrease of t^J as the decimation process pro- 
ceeds. After a large number of decimation steps, the Lya- 
punov exponent is asymptotically related to the effective 
hopping tg 6 ^/" 1 in the following way: 



N 



Aeff) 
L 0,N 
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(N-2), m 
e N-l \E) 



(E) = t 



(e//) 
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N-1,N 



E-e$- 2) (E) 



,(6) 



(7) 



{eff \E)\}. 



III. RESULTS AND DISCUSSION 



where, after N — 1 decimations, e 



N-l 
N 



denotes the renor- 

(eff) 



malized diagonal element at site TV and t Q ^ indicates 
the effective renormalized hopping connecting the sites 
and N . The behavior of the effective interaction t^ff^ 
during the decimation procedure can be used as another 



In the present study, we use chains with 10 7 sites in 
order to calculate the Lyapunov exponent within the 
entire energy band. We start our analysis investigat- 
ing the main features of a diluted ID Anderson model 
with segments consisting of two sites with identical on- 
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FIG. 2: a) Lyapunov exponent 7 versus E for a chain with 
N = 10 7 sites diluted by dimers with eo = 1.0 . We found that 
7 is finite for all energies except at the two resonance energies 
where 7 vanishes linearly (see b) and c)). For eo = 0.0, one 



E, 



\(2/3) 



has a slower non-linear vanishing of 7 oc {E 
long-dashed line in fig b). (d) The effective hopping t^f/- 1 
displays an oscillating behavior as a function of N for E = 
0.0 (extended state) and an exponentially decaying behavior 
for E — 0.1 (localized state). 



site energies eo- The emergence of extended states for 
this particular dilution was analytically demonstrated in 
Ref [H i3. The density of states (DOS) for general 
ID tight-binding models can be obtained by using the 
negative eigenvalue theorem[35l|. In Fig. 1(a) we show 
typical plots of the DOS of chains with N = 10 5 sites 
and several eo values. In this diluted chain, the density 
of states displays two pseudo-gaps for eo > 0. The widths 
of these pseudo-gaps increase as a function of eo- These 
pseudo-gaps are reminiscent of the gaps present in the 
absence of disorder and their width is proportional to eo 
for any length of the diluting segment (see fig. lb) , thus 
following: a similar trend of the gap found in the two band 
model 37]. The presence of disorder rounds one of the 
band edges which displays an exponentially decaying tail. 
On the other hand, the opposite band edge remains as a 
singularity of the DOS. The surviving of this singularity 
in the presence of disorder will be shown to influence the 
wave-packet dynamics of the diluted Anderson model. 

In Fig. 2(a) we depict the Lyapunov exponent 7 versus 
E for chains with N — 10 7 sites and eo = 1.0. We calcu- 
late 7 using both transfer matrix and decimation renor- 
malization group approaches. The results were identical 
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FIG. 3: a) The Lyapunov exponent 7 versus E for an An- 
derson chain diluted by segments with k = 3 sites exhibiting 
specular symmetry f(l) = {£2,^3, £2} where €2 = 0.8 and 
£3 = 1.0. All calculations were performed for a chain with 
N = 10 7 sites, (b) The effective hopping tg(/ as a function 
of N for the resonance energies E c = {0.8, 2.31774, -0.51774} 
corresponding to extended states. 



within our numerical accuracy. We found that j(E) is 
finite for all energies except at the two resonance energies 
(E = 0.0 and E = 2.0) where the Lyapunov exponent is 
of the order -A.. Therefore, all states are exponentially 
localized except at these two resonance energies where 
7 vanishes in the thermodynamic limit. This result co- 
incides with the analytical prediction of extended states 
at E = e ± 1 .24]. In Fig. 2(b) and 2(c), we point out 
the linear vanishing of 7 around E = and E = 2.0 (see 
full lines in fig Fig. 2(b) and 2(c)). For e = 0.0, we 
found a slower non-linear vanishing of 7 oc (E — i? c )( 2 / 3 ) 
(see dashed line in fig. 2(b)). in perfect agreement with 
ref. [l^. The larger exponent found for eo ^ can be 
attributed to fact that the resonance energies for eo ^ 
are precisely at the DOS band edge singularities, which 
are absent for eo = 0. To further characterize the ex- 
tended nature of these resonant states, we plot in Fig. 2d 



the effective interaction t 



(eff) 
0,N 



versus N. For the resonant 



state at E — 0.0 the Lyapunov exponent 7 vanishes due 
to the oscillating behavior of t^ffi . The extended char- 
acter of this state is, therefore, reflected by a finite effec- 
tive hopping amplitude between the sites located at the 
chain ends. For E = 0.1 the effective hopping decreases 
exponentially, indicating the localized character of non- 
resonant states. We have also considered longer diluting 
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FIG. 4: (a) The Lyapunov exponent 7 versus E for an 
Anderson chain with N = 10 7 sites diluted by segments 
with k = 3 sites exhibiting the distance to center symmetry 
f(l) = £; = {0.8,0.7,0.6}. In this case, there is a single ex- 
tended state located at the resonance energy E = 0.7. (b) The 
effective hopping t^ff as a function of N for E = 0.7 shows 
the oscillating behavior characteristic of extended states. 



sequences with constant on-site energy. The number of 
extended states always corres pon ds to the number of sites 
in the sequence, as expected [231 l24j . 

From the theoretical point of view, the basic condition 
to find delocalized states in diluted disordered systems is 
the existence of some symmetries in the periodic function 
f(l) which defines the site energies e; = /(/) within the 
diluting segment. The relevant symmetries are|24j: 

a) Specular reflection symmetry with respect to center of 
the diluting segment, namely, f(l) — f(k + 1 — 1), with 
I = 1,2, ...k. 

b) Distance to the center symmetry, namely, \f(lo) ~ 

= \f(h)-f(k+j)\, with j ranging from j = 1 
to j = (k — 1) /2 and lo being the position of the central 
segment site. This kind of symmetry can be present only 
in segments with an odd number of sites. 

Delocalized states emerge whenever the diluting func- 
tion exhibits one of the above symmetries. In the follow- 
ing, we consider two examples using diluting segments 
with either one of the above symmetries present. 

To observe the emergence of extended states in the 
diluted Anderson model with specular reflection symme- 
try, we considered a diluting segment with k = 3 sites: 
f(l) = {£2, £3, £2}, where the site energies are £2 = 0.8 
and £3 = 1.0, respectively. In fig. 3(a) we show the Lya- 
punov exponent 7 versus E. All calculations were per- 



formed for N = 10 7 sites. After executing the renormal- 
ization process, we found that the Lyapunov exponent 
vanishes for E c w {0.8, 2.31774, -0.51774} , in agreement 
with the block decimation result [24|. In addition, we 
show in Fig. 3(b) the oscillating pattern exhibited by the 

effective hopping t^J^ as a function of N for the same 
critical energies, thus confirming the extended nature of 
these states. 

A similar analysis can be made to investigate a diluted 
Anderson chain with the distance to the center symme- 
try. We considered, for this case, the following diluting 
function: f(l) = ej = {0.8, 0.7, 0.6}. From the theoreti- 
cal calculation of [24|, one can anticipate that this chain 
will have only one extended state with energy E c = 0.7 . 
Fig. 4(a) and Fig. 4(b) clearly show that the Lyapunov 
exponent vanishes only for E c = 0.7 and the effective 
hopping t^J^ displays an oscillating behavior at this 
critical energy. 
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FIG. 5: The scaled mean-square displacement a 2 (t)/t ' 5 ver- 
sus times t for the ID diluted Anderson model with 7 x 10 4 
sites for: a) Constant diluting function eo = f(h) = eo = 1 
with l\ = 2, 20. (b) Diluting function with specular reflec- 
tion symmetry f(l) = {€2,63,^2} with £2 = 0.8 and £3 = 1.0 
(c) Diluting function with the distance to center symmetry 
f(l) = {0.8, 0.7, 0.6}. A sub-diffusive behavior a 2 oc t ' 5 takes 
place for any symmetry supporting extended states. 



In order to investigate the influence of the above res- 
onant extended states on the transport properties of the 
diluted Anderson model, we calculate the time depen- 
dence of the mean-square displacement of an initially 
localized wave-packet. Starting with one electron fully 
localized at the Anderson site closer to the chain center, 
the wave-function amplitudes b n (t) were obtained from 
the following equations of motion: 

ib n {t) = e n b n (t) + (b n -i(t) + b n+1 {t)) n=l,2,...,N. 

(8) 
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FIG. 6: Dashed lines represent the normalized density of 
states (DOS) for the diluted chain with a) specular reflection 
symmetry f(l) = {0.8, 1.0,0.8} and b) distance to the center 
symmetry f(l) = {0.8,0.7,0.6}. In both cases, the resonance 
energies associated with extended states are coincident with 
DOS singularities. The DOS displays pseudo-gaps for eo > 0, 
which are reminiscent of the gaps present in the DOS of the 
corresponding pure models (shown as solid lines). 

Using a Runge-Kutta algorithm, we solve the above set 
of coupled equations and calculate the mean-square dis- 
placement 

* 2 (t) = J2(n-no) 2 \b n (t)\ 2 . (9) 

n 

To minimize end effect, our numerical calculation was 
performed in a very large chain with N — 7 x 10 4 sites. 
As a consequence, the site amplitudes at the ends of the 
chain (&i(t) and 6jv(i)) are always negligible. The numer- 
ical integration of equations (8) was performed using the 
Runge-Kutta algorithm with precision At smaller than 
1CP 2 . In Fig. 5(a-c), we exhibit our main results for the 
time evolution of the mean-square displacement cr 2 (t) on 
diluted chains considering all relevant symmetries of the 
diluting function. In all cases, we obtained a clear sub- 
diffusive behavior: a 2 ot t 0,5 . This sub-diffusive spread of 
the wave-packet is related to the fact that the resonance 
energies supporting extended states occur at band edge 
singularities of the DOS, as depicted in Fig. 6. A simi- 
lar diffusive-like spread of the wave-packet has also been 
observed to occur with collective excitations in other dis- 



ordered systems with the Lyapunov exponent vanishing 
at DOS singularities as, for example, in random harmonic 
chains [34| and disordered ferromagnetic chains 0, ■ 



IV. SUMMARY AND CONCLUSIONS 

In summary, we investigated the ID diluted Anderson 
model where every Anderson impurity is diluted by a set 
of site energies given by a diluting function f(l). Us- 
ing the renormalization group approach and the trans- 
fer matrix technique, we obtained the Lyapunov expo- 
nent 7 in several diluting cases exploring distinct sym- 
metries supporting extended states, such as constant di- 
luting function, diluting function with either specular 
reflection or distance to the center symmetry. In all 
cases, we identified the resonance energies at which the 
model exhibit extended states, in full agreement with 
previous results [2^, |24|. We also studied the tempo- 
ral spread of an initially localized electron wave-function 
in these diluted chains by following the time dependence 
of the wave-packet mean-square displacement. We found 
that, associated with the presence of extended states at 
resonance energies located at DOS band edge singular- 
ities, the electron wave-packet displays a sub-diffusive 
spread(er 2 (i) oc i 5 ). This result is consistent with the 
general picture that the wave-packet dynamics depends 
on the relative location of resonance energies and DOS 
singularities [36| . 
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